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Abstract 

We have observed coherent oscillations in a multi-level quantum system, formed by a current- 
biased dc SQUID. These oscillations have been induced by applying resonant microwave pulses of 
flux. Quantum measurement is performed by a nanosecond flux pulse that projects the final state 
onto one of two different voltage states of the dc SQUID, which can be read out. The number of 
quantum states involved in the coherent oscillations increases with increasing microwave power. 
The dependence of the oscillation frequency on microwave power deviates strongly from the linear 
regime expected for a two-level system and can be very well explained by a theoretical model taking 
into account the anharmonicity of the multi-level system. 

PACS numbers: 03.67.Lx,05.45.-a,85.25.-j 
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Up to now, in view of quantum information processing, experiments in solid state devices 
have concentrated only on the implementation of two-level quantum systems. A variety of 
quantum circuits based on Josephson junctions □ □DO and quantum dots |5( have been 
proposed and realized. Rabi oscillations have been observed showing that the two lowest 
levels can be manipulated coherently. However, the two-level system is not the only one 
useful for quantum computation. For instance, multi-level systems are of interest for solving 



database search problems using quantum algorithms (J, as was demonstrated with Rydberg 
atoms [7|. Recent theoretical proposals discuss the use of multilevel quantum systems in solid 



state devices for quantum information processing In superconducting devices such as a 
current-biased Josephson junction [9] or a rf SQUID Q], many experiments were performed 
demonstrating the multi-level quantum nature of these systems. However no experiment 
has probed coherent behaviour yet. In this letter, we report, to our knowledge, the first 
observation of coherent oscillations in a multi- level solid-state circuit, based on Josephson 
junctions. The non-linearity of the Josephson junction plays a crucial role in the quantum 
dynamics of the device. 

Specifically, the quantum system that we study in our experiments is a current-biased dc 
SQUID. It consists of two Josephson junctions each with critical current Iq and capacitance 
Co. The junctions are embedded in a superconducting loop of inductance L s , threaded by a 
flux $. Since L S I « $ /27r, the phase dynamics of the two junctions can be treated as that 
of a fictitious particle of mass m = 2C , n($n/27r) 2 moving in a two-dimensional potential, 
which contains valleys and mountains (jjj, llJ, Here $o — h/2e is the superconducting 



flux quantum. The local minima are separated from each other by saddle points, where 
the particle can escape. Along the escape direction, the potential is cubic and can be 
characterized by a frequency uj p and a barrier height AU, see Fig. [TJi. These two quantities 
depend on the magnetic flux and vanish at the critical current I c of the SQUID. We assume 
a complete separation of the variables along the escape direction and the transverse one 
by neglecting the coupling terms between these two directions. In this approximation, the 
dynamics of the SQUID's phase along the escape direction is similar to the dynamics of the 
phase of a current-biased single Josephson junction jl^l 15^ . The parameters I c , AU and u p 
are renormalized, thereby taking into account the two-dimensional nature of the potential, 
as it was demonstrated in jll|, [l^. For bias currents < J c , the particle is trapped in a 
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local minimum and its quantum dynamics is described by 

H = hiuj p {P 2 + X 2 ) - ahu p X 3 , (1) 

where P = (1/ ^mftw p )P and X = {^Jmuo p /K)(j) are the reduced momentum and position 
operators, respectively. Here, P is the operator conjugate to (ft; a = l/(6a)[2(l — ib/%)] ~ 5 ^ 8 
is the relative magnitude of the cubic term compared to the quadratic (harmonic) term. The 
parameter a is a constant, a ~ 11 for our SQUID. For lb well below I c , AU ^> hu p and many 
low-lying quantum states are found near the local minimum. These states, describing the 
oscillatory motion within the anharmonic (cubic) potential, are denoted \n) for the nth level, 

with n — 0, 1, 2, The corresponding energies E n were calculated in M .Q.«i ng 

through the barrier can be neglected for the lowest-lying states with nhu p <C AU. The effect 
of a time-dependent flux $(t) can be included in Eq. ((TJ) by adding the time-dependent term 
a(t)hu p X where a(t)is proportional to $(t) [16]. In our experiment, deep quantum states 
are excited by applying microwave (MW) flux pulses characterized by their frequency u, 
amplitude &mw an d duration AT. Starting with 1^(0)) = |0), the state evolves with the 
MW pulse of duration AT into a coherent superposition |\1>(AT)) = a |0) + ai 1 1) + . . . + 
a n \n) + ... 

Our procedure to perform quantum experiments consists of four successive steps as de- 
picted in Fig. ^p. A bias current h is switched on through the SQUID at fixed magnetic flux 
$d c to prepare the circuit at the working point in the initial state |0). The application of a 
MW flux pulse produces the superposition |\1/(AT)). Then, a flux pulse of nanosecond dura- 
tion is applied which brings the system to the measuring point such that the state |^(AT)) 
is projected onto either the zero or the finite voltage state of the SQUID (Fig. \Bp). The 
result of this quantum measurement can be read out by monitoring the voltage V ont across 
the dc SQUID. Finally, h is switched off such that the circuit is reset and the experiment 
can be repeated to obtain state occupancy. 

The aforementioned quantum measurement procedure was discussed theoretically by us 
m Ref. 0- Differently from previous experiments which used current bias pulses |2j, y] 
or MW pulses we implemented a quantum measurement using a large-amplitude flux 
pulse with nanosecond duration. This flux pulse reduces the SQUID critical current to a 
value very close to the bias current such that AU ~ hjj p . Ideally, the pulse with optimal 
amplitude projects the excited states with n > onto the voltage state (V ont is twice the 
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superconducting gap); the ground state |0) is projected onto the zero- voltage state. As the 
SQUID is hysteretic, the zero and finite volta ge s tates are stable. The efficiency of this 
one-shot measurement is estimated to be 96% [16j. In the experiment, we use a pulse of 
2.5ns duration whose 1.5ns rise and fall times are long enough to guarantee adiabaticity: 
transitions between |0) and excited states are estimated to occur with a probability less than 
1% for a typical pulse amplitude of 0.06$o- The measurement pulse can be applied with a 
variable delay after the end of the MW pulse. In the measurement procedure, this delay is 
kept as short as 1.5ns to limit relaxation processes. The result of the quantum measurement 
is obtained by measuring the escape probability P esc repeating the experiment up to about 
4000 times. 

The measured dc SQUID consists of two aluminum Josephson junctions with I = 
3.028/M and Co = 0.76pF, coupled by an inductance L s = 98pH (see Fig. |T]i). The quan- 
tum circuit is decoupled from the external classical circuit by long on-chip superconducting 
thin wires of large total inductance, L e = 15nH, and a parallel capacitor, C e = 150pF. The 
quantum circuit and the long superconducting wires are realized by e-beam lithography and 
shadow evaporation. The nominal room temperature microwave signal is guided by 50Q 
coax lines and attenuated twice by 20dB (at 1.5K and 30mK, respectively) before reaching 
the SQUID through a mutual inductance M s = 0.7pH. The MW line is terminated by an 
inductance estimated to be Lmw — InH. Special care was taken to avoid spurious envi- 
ronmental microwave resonances. The chip is mounted in a shielded copper cavity, cooled 
down to about 30 mK, whose cut-off frequency is above 20 GHz. Moreover, the thin film 
capacitor, C e , close to the chip decouples the microwave signal from the low-frequency bias 
lines. 

The escape probability is first obtained by measurements using long (duration At = 
50/is) pulses of the bias current From the dependence of the escape current (defined as 
the current at which P esc = 1/2) on $^c, the experimental parameters of the SQUID are 
extracted llj|. At the maximum value of I c obtained for / $n = —0.085, the measured 
escape rate can be fit by the well-known MQT formula Q| yielding the same SQUID 
parameters within a 2 % error. This result confirms that the circuit remains in the ground 
state when only I b pulses are applied. For other values of $d c , the width of the escape 
probability is larger than the MQT prediction by a factor up to 4 indicating a residual 
low- frequency flux noise in our sample This noise reduces the efficiency of our quantum 
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measurement since it smears the escape probability difference between the ground state and 
the excited states. In the experiment the nanopulse amplitude is adjusted in order that 
the escape from |0) is close to 1%. For such a pulse amplitude at 3>dc/®o = 0.095, we 
have estimated that the escape probabilities of the excited states |1), |2), |3) and higher 
states are about 30 %, 60 %, 90 % and 100% respectively. The sensitivity to flux noise 
is weakest at the "optimal" point $dc/$o — —0.085; however, the SQUID's sensitivity to 
MW and measurement pulses is also weak such that this point cannot be used for accurate 
measurements. 

Spectroscopic measurements were performed by sweeping the frequency of a MW flux 
pulse of 25ns duration in the low power limit. In the inset of Fig. El we show the cor- 
responding resonance peak found for P esc versus microwave frequency v. The resonant 
frequency z/oi depends on If, as shown in Fig. El this dependence can be very well fit by 
the semiclassical formula for the cubic potential [15]. The parameters extracted from this 
fit are consistent within 2% error with the parameters extracted from the critical current 
versus flux dependence or from the MQT measurements at $d c /$o — —0.085. The large 
width Auqi = 180MHz of the resonance peak is consistent with the presence of a residual 
low-frequency flux noise. Indeed, since the flux through the SQUID fluctuates slowly, the 
expected resonant frequency changes from one measurement to the other. 

We also measured the lifetime of the first excited states, analyzing the decay of the 
resonance height. Upon increasing the delay time between the end of the MW pulse and 
the nanosecond dc measurement pulse, the peak height is found to decay with times ranging 
from 14ns to 60ns. The frequency broadening associated with these lifetimes is estimated 
to be smaller than 10MHz, i.e., not only less than the width Az/ i of the resonance peak but 
also less than the detuning frequency associated with the anharmonicity between adjacent 
levels. 

By applying short MW pulses at the resonant frequency v = uqi, we induce coher- 
ent quantum dynamics. Rabi-like coherent oscillations were observed by measuring escape 
probability versus MW pulse duration AT. In Fig. the escape probability oscillates at 
a frequency Q co h/2n of about 300MHz. This frequency increases as the MW flux power 
increases and ranges from about 100 MHz to 1000 MHz in our experiment. It is always 
much smaller than the resonance frequency v m . Oscillation amplitudes as large as 70% were 
observed for the largest MW power. The oscillations are damped with a characteristic at- 
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tenuation time of about 12ns. Similar coherent oscillations have been observed at different 
working points. The dependence of fl co h as a function of the MW amplitude $mw is shown 
in Fig. |U The linear dependence predicted by Rabi theory for a two level system [l^ is not 
observed in our measurements. 

To analyse the observed coherent oscillations, we use the following model. We ignore 
relaxation and decoherence processes and consider the applied MW frequency to be the 
resonant frequency z^oi- If we furthermore assume the MW pulse to couple the eigenstate 
|n) to its nearest-neighbor levels only, the time-dependent part of the Hamiltonian reads 
a(t)huj p X « 6^^cos(2vrz/ it)^p^ n ^/n/2[\n)(n - 1| + \n - l)(n|], where b ~ 34 at the 
considered working point. As £l co f l /2ir <C uqi, we can use the rotating wave approximation 
and treat this term non-perturbatively. We obtain the coefficients a n of the coherent super- 
position |^(AT)) generated by a pulse of duration AT: a n = J2 k (ek\0) (n\ek)e~ tXkAT ^ h , where 
\ek) and Afc are the eigenstates and eigenvalues of the Hamiltonian in the rotating frame. The 
latter are determined numerically and then used to calculate the probabilities \a n (AT)\ 2 , 
whose oscillatory time dependence is determined by the frequency differences — Xj. At 
very low MW amplitude (b^Mw/^o^p <^ 2£?i — (E 2 + Eq)), we obtain the two-level limit: 
|^(AT)) oscillates between |0) and |1) at the Rabi frequency flcoh — bu p /\/2^ MW /^Q jl7|. 
At larger MW amplitude, an increasing number of states is predicted to participate in the 
oscillations Q. As an examp i e! in Fi g. a the theoretical escape probability is plotted as 
function of MW pulse duration taking into account the finite efficiency of our measurement. 
The MW amplitude &mw/®o — 0.002 was adjusted to fit the measured oscillation frequency. 
The MW power calibration found from the fit is consistent, within a 30 % error, with the 
actual applied MW amplitude and with the microwave line and coupling parameters, which 
were measured independently. The model calculation predicts very large amplitudes, con- 
sistent with the experimental results. At AT = 1.6ns the amplitude is close to 100%. In 
Fig- Et we indicate the occupancies |a n | 2 of the states |n) that participate in the coherent 
superposition state at this time. We see that it is concentrated mainly on the states |4), |5) 
and 1 6). Finally, the model predicts a slight beating in the oscillations at AT ~ 20ns where 
the amplitude reaches a minimum. These beating phenomena will not be analyzed in our 
experiment because they occur at time scales of the order of the above-mentioned attenu- 
ation time. For short duration times, the predicted coherent oscillation frequency is given 
by min fcj j{(Afc — \i)/h}. In this limit and using the MW calibration deduced from Fig. |HJ we 
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calculate O co ft,/27r as a function of $mw/®q- The result is plotted in FigJU together with the 
experimental one. The perfect agreement between them shows that our model captures the 
physics of the coherent oscillations. 

We reported, to the best of our knowledge, the first observation of coherent oscillations 
in a multi-level solid-state-based integrated circuit. In our experiment, we have induced 
coherent superpositions of quantum states which involves many levels using monochromatic 
microwaves. In order to perform the quantum measurement, we implemented a new mea- 
surement procedure based on nanosecond flux pulse through the SQUID. The agreement 
between our results and the theory proves that the essential physics of our non-linear quan- 
tum circuit is well understood. 
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FIG. 1: Quantum experiment and measurement procedure in a dc SQUID, (a) Illustration of 
the fictitious particle trapped in a deep anharmonic well at the working point, (b) Experimental 
procedure, as explained in the text, (c) The well at the measuring point: during the nanosecond 
flux pulse, the particle can escape or remain in the well depending on its quantum state, (d) 
The dc SQUID consists of two identical Josephson junctions coupled by an inductance L s . The 
quantum circuit is decoupled from the environment, symbolized by Z(u), by a large inductance 
L e and a capacitor C e ; it is coupled by a mutual inductance M s to the MW pulse or nanosecond 
pulse through a 5012 coax line, terminated by an inductance Lmw- 
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^ 4.4 4.6 4.8 5 

Current bias I b (|uA) 

FIG. 2: Resonant transition frequency foi between states |0) and |1) as function of the bias current 
at $rf c /$o = 0.095. The experimental dependence of vq\ (symbols) on is perfectly described 
by the semiclassical model |l5( (solid line). Spectroscopy is performed by measuring the escape 
probability P esc versus applied frequency v (points in inset) at microwave power A=-48dBm where 
A is the nominal microwave power at room temperature. The fit to a Gaussian (solid line) defines 
the resonant frequency v§\ = HA3GHz and the width Ai/qi = 180MHz. 
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Microwave pulse duration (ns) Level index, n 

FIG. 3: a) Measured escape probability versus MW pulse duration at the working point, denned 
by &dc/&o = 0.095 and I\, = 4.7/jA, with about 15 levels localized in the anharmonic well, b) 
Predicted escape probability for a MW amplitude &mw/&o = 0.002. c) a n coefficients of the state 
obtained after 1.6ns duration of the MW pulse corresponding to the first maximum of the coherent 
oscillations. 
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Microwave amplitude $mw /^o 



FIG. 4: Frequency of coherent oscillations versus the normalized MW flux amplitude. Dots are 
experimental results, the solid line the theoretical prediction for a multi-level anharmonic quantum 
system. Dashed line is the Rabi theory for a two-level system. Experimental calibration of &mw 
is deduced from the fit and is consistent with the estimated %if applied to the set-up. 
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